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We study the properties of one-dimensional gapped Heisen- 
berg antiferromagnets in the presence of an arbitrary strong 
staggered magnetic field. For these systems we predict a uni- 
versal form for the staggered magnetization curve. This func- 
tion, as well as the effect the staggered field has on the energy 
gaps in longitudinal and transversal excitation spectra, are 
determined from the universal form of the effective potential 
in 0(3)-symmetric 1+1-dimensional field theory. Our theo- 
retical findings are in excellent agreement with recent neu- 
tron scattering data on ifeBaNiOs (R = magnetic rare earth) 
linear-chain mixed spin antiferromagnets. 



One-dimensional isotropic Heisenberg antiferromag- 
nets with an exchange gap in the magnetic excitation 
spectrum have been at the center of theoretical and ex- 
perimental attention for almost two decades. This class 
of materials includes integer-spin Heisenberg chains |l]] 
(commonly referred to as Haldane-gap systems) and half- 
integer spin ladders with an even number of legs Q . Due 
to the presence of strong quantum fluctuations in these 
systems the staggered magnetization has a finite correla- 
tion length. The principal feature of the excitation spec- 
trum is a degenerate triplet of sharp spin-1 excitations 
commonly referred to as magnons, separated from the 
ground state by a finite energy gap A. In recent years 
much work was aimed at understanding the behavior of 
such gapped ID antiferromagnets in the presence of an 
external uniform magnetic field However, the ef- 

fect of a staggered field, that couples directly to the order 
parameter of the classical system, has not been inves- 
tigated in sufficient detail. This is mainly due to the 
fact that a strong magnetic field modulated on the mi- 
croscopic scale was thought to be all but impossible to 
realize experimentally ||. A breakthrough came with 
neutron scattering experiments on /i^BaNiOs (R = mag- 
netic rare-earth) linear-chain nickelates and their inter- 
pretation in terms of non-interacting Haldane spin chains 
immersed in a strong effective staggered exchange field 
(@,|). In _R 2 BaNi05 compounds almost perfectly isotropic 
antiferromagnetic 5=1 chains are formed by the Ni 2+ 
ions. The effective staggered field is generated by the 
R 3+ sublattice that becomes ordered magnetically below 
some Neel temperature T)v. The staggered field intensity 
is proportional to the magnitude of the ordered moment 
of i? 3+ magnetic ions and can be controlled in an experi- 



ment indirectly, by varying the temperature. One of the 
most significant results was the first direct measurement 
of the staggered magnetization curve M S (H S ) of a Hal- 
dane spin system |)j . It was found that of particular value 
as model systems are (NdzYi-^BaNiOs species, where 
the effective interaction between R and Ni-sublattices is 
of Ising type jlO| , so that the transverse excitations on 
the Haldane chains are effectively decoupled from those 
on the R- subsystem. Thus the effect of the staggered 
field on transversal magnetic excitation spectrum of an 
isolated Ni 2+ Haldane chain could be measured experi- 
mentally @. 

In H we gave a qualitative theoretical description of 
quantum disordered antiferromagnets in the presence of 
a staggered magnetic field and discussed the results in 
relevance to existing data on ^BaNiOs materials. The 
principal conclusion was that in a weak staggered field 
the energy gap A increases in proportion to the square 
of induced staggered moment on the Haldane chains. It 
was also shown that a staggered field partially lifts the 
degeneracy of the magnon triplet, the gap in the longitu- 
dinal mode being three times more sensitive to H s than 
that in two transversal magnons. In the present paper we 
refine this approach and obtain quantitative predictions 
that we directly compare to recent experimental data for 
(Nd x Yi_ x )2BaNi05. We demonstrate that the behav- 
ior of both transversal and longitudinal energy gaps in 
the presence of a staggered field is contained in the stag- 
gered magnetization curve M S (H S ). Our central result 
is that for a variety of gapped one-dimensional antiferro- 
magnets M S (H S ) has a universal shape defined by three 
experimentally accessible parameters: zero-field magnon 
energy gap A, spin- wave velocity v, and the renormal- 
ization constant Z, related to the residue of the magnon 
pole at H s = 0. 

A traditional theoretical description of spin dynamics 
in one-dimensional quantum antiferromagnets is based 
on the mapping of these systems to the 1+1-dimensional 
0(3) Non-Linear Sigma Model (NLSM) with one spa- 
tial and one temporal coordinates. This approach was 
first introduced by Haldane Q for a Heisenberg antifer- 
romagnetic spin chain, and later extended to a variety 
of other systems (for a recent review see [0). Although 
this mapping is based on a quasiclassical approximation, 
and in principle it should work well only for large spins 
S 3> 1, it gives quantitatively correct predictions for any 
spin value. In the absence of topological term the NLSM 
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Lagrangian in the presence of external (staggered) field 
can be written as 
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Here is a three component unit vector (0 2 = 1), point- 
ing in the direction of the local staggered moment, vq 
is the bare spin wave velocity of the system, and go 
is the dimensionless parameter, controlling the strength 
of quantum fluctuations. Notice that in our notation 
the staggered field H s is coupled to the local staggered 
spin. The staggered magnetic field Hs m \ coupled to lo- 
cal staggered magnetic moment, is proportional to H s : 
g m fJ,BHs — H s , where g m is the g-factor of a magnetic 
ion, and \xb is the Bohr magneton. In a large-S map- 
ping of the Heisenberg spin chain to the NLSM yj] one 
has vq = 2JS, and go = 2/5. In a more general situa- 
tion these parameters should be treated as phenomeno- 
logical constants fine-tuned to give the correct low en- 
ergy properties of the system. The 1+1-dimensional 
0(3) NLSM is always in a disordered state. The cor- 
relation length in the absence of staggered field is given 
by £ ~ a exp(27r/<7o), where a is the lattice spacing. The 
gap A in the excitation spectrum is related to the corre- 
lation length via a usual relation A = hv/£. 

In order to calculate the macroscopic properties of 
the system such as its staggered magnetization curve 
or the excitation spectrum near the AFM zone cen- 
ter, one needs to coarse-grain the NLSM Lagrangian 
by integrating out the large q and u> (small x and t) 
degrees of freedom. The change of parameters of the 
Lagrangian as a result of this coarse-graining is de- 
scribed by the Renormalization Group (RG) flow equa- 
tions. In real space the renormalization procedure cor- 
responds to replacing the field 0(x,t) with r (x,t) = 

( v / 12 ) I\ X '-x\<i,\t'-t\<i/v'?( x '' t ') dx ' dt '- lt is eas y to see 

that the coarse-grained field variable r no longer has a 
well defined length. In other words, as a result of coarse- 
graining the NLSM Lagrangian, characterized by a rigid 
constraint 1 — 1, becomes a "soft-spin" Lagrangian. 
The length of coarse-grained field variable r has a prob- 
ability distribution, defined by some effective potential. 
The fully renormalized Lagrangian can be written as 



£ = dx 



' hit) -lit] -^-^ 



(2) 

In this expression we have introduced the renormalized 
staggered field variable 0, defined through \J~Z§ = S0. 
The renormalization parameter Z was fine-tuned to give 
the desired form of the derivative terms, with v being the 



true spin wave velocity. To specify this fully renormal- 
ized Lagrangian we need to know the form of effective 
potential C/(|0|). For small (j> it can be expanded in Tay- 
lor series in <p 2 . The quadratic term can be written as 
(A 2 /2v)cf> 2 , where A is the true energy gap of the magnon 
spectrum. Indeed, the quadratic part of (||) should de- 
scribe a triplet of "relativistic" non- interacting bosons. 
The spin wave velocity v plays the role of the speed of 
light in its relativistic analogue. The above form of the 
quadratic term correctly reproduces the spin-wave dis- 
persion £(ir + k) = VA 2 + v 2 k 2 . The relativistic analogy 
also determines the intensity of the magnon pole in the 
spin correlator as 
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In 1+1 dimensions it is convenient to write the Taylor 
expansion of ?7(|<^|) in terms of the set of dimensionless 
parameters um defined by 

U{\$\) = ^{\\$\ 2 + + ^u 6 \$\ e + ■■■) (4) 

The effective potential, truncated at the \4>\ 4 term, corre- 
sponds to the Ginzburg-Landau Lagrangian introduced 
by Affleck [[l2| on phenomenological grounds to describe 
Haldane-gap systems. 

The linear staggered susceptibility at zero external 
staggered field can be derived from the quadratic part of 
Lagrangian describing non-interacting magnons. Indeed, 
the expectation value of 4> in the presence of a weak stag- 
gered field is easily obtained by balancing the quadratic 
term and the source term in Eq.(^). The result is the 
single mode contribution to zero-field staggered suscep- 
tibility 
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We can compare this prediction to the numerical results 
for the S=l Heisenberg chain. Using the numerical values 
Z = 1.26 (g in their notation), v — 2.49J, and A = 
0.41J, reported in |^3|, the single mode approximation 
gives x^(0) — 18. 7/ J in excellent agreement with Monte 
Carlo result 21(1)/ J [[u) . This is a manifestation of the 
well known fact that in Haldane-gap systems virtually all 
spectral weight at the AFM zone center is concentrated 
in the magnon triplet. 

The main concern of the present paper is the non- 
linear behavior in arbitrary strong staggered fields. In 
order to describe these effects quantitatively one needs 
to know the numerical values of dimensionless couplings 
U2n, which in principle should depend on the parame- 
ters of the bare Lagrangian. If, however, the correla- 
tion length is sufficiently long one can safely assume that 
these couplings are at their RG fixed point values. For- 
tunately, very accurate numerical values of the universal 
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fixed point couplings 1*4, Uq, and were recently ob- 
tained by Pelissetto and Vicari in |l5|,[l6| . They carefully 
compared the results of e, l/N, high temperature, and 
strong coupling expansions of 0(A^)-symmetric models in 
d dimensions with the results of Monte Carlo simulations. 
For d — 2, N = 3 their estimates are 114 — 11.8(1), u@ — 
3.33(10) xul = 460(20), u 8 = 20(5) xv% = 33,000(8,000) 
p7[ . In pl| the deviations of couplings from their fixed 
point values were estimated as «2n(£) — U2n(oo) ~ l/£ 2 - 
For the S — I Heisenberg antifcrromagnet the correlation 
length £ ~ 6 (in units of lattice spacing). Therefore, one 
could reasonably expect the deviation of u% n from their 
fixed point values to be around 1/36 ~ 3%! This fact is 
confirmed numerically in Fig. 2 of fl5}| . 

From the above it follows that the quartic term A|</>| 4 
used by Affleck Q to describe the effects of pairwise 
magnon repulsion has in fact a universal strength A = 
(u 4 /4!)(A 2 /i;) ~ 0.49A 2 /V In their theoretical study 
of the effect of the external field on the excitation spec- 
trum in NENP Mitra and Halperin made a rough 
estimate of the value of A. By matching the first term in 
perturbative l/H s high field expansion of the magneti- 
zation function with a small H s expansion they have got 
an order of magnitude estimate, which in our notation 
corresponds to A ~ ZA 2 /Av ~ 0.3lA 2 /t> rather close to 
our more refined result. 

The effective potential manifests itself in the staggered 
magnetization curve of a ID gapped antiferromagnet. 
Let us select z-axis along the external staggered field. 
The expectation value of the field (f> is defined by the 
minimum of the total potential energy (effective potential 
plus an external field term) located at \fZH s = U'(((f> z )). 
Therefore, the staggered magnetization curve M S (H S ) is 
defined by the equation 
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Using the numerical estimates for 114, uq, u% quoted 
above, and Z = 1.26 from [|lj|, for S — 1 chain we ob- 
tain x {s) (0)H s = M s (l + 1.56M 2 + 2.4A/ 4 + 3.27M S 6 ) 
This result can be directly compared to the M S (H S ) 
curve measured experimentally in NdzY^^BaNiOs || 
(Fig. 1). The solid line is a fit to the experimental data 
with x^(0) being the only adjustable parameter. We 
note that the result of the fit - Xcxp (0) = O.SSmeV- 1 
differs from the expected value Xtheor(^) = 18. 7/ J ~ 
0.85meV -1 . It has to be emphasized however, that the 
experimental scaling of the abscissa in Fig. 1 in units of 
magnetic field heavily relies on a series of assumptions 
and simplifications regarding the properties of rare earth 
ions, and the exact applicability of the mean-field model 
in a wide temperature range . The experimentally de- 
termined shape of M S (H S ), on the other hand, is a more 



robust result and is in excellent agreement with our the- 
oretical predictions. 

The effective potential also contains information on 
the behavior of the excitation spectrum in the pres- 
ence of a finite staggered field. The magnon exci- 
tations are by definition deviations of local staggered 
magnetization from its equilibrium value. When the 
0(3) symmetry is broken by the external field, the de- 
generacy of the magnon triplet is partially lifted. In 
this case one should distinguish between the gap A|| 
in the longitudinal branch of 5<j) z excitations and the 
transversal gap Aj_ in the doublet of 5<fi x , 8<f> y excita- 
tions. Both these gaps are determined by the quadratic 
terms in the expansion of U(\(j) + 6<fr\) ~ £^(H) + 
^(A 2 /2w) (5</> 2 . After some straightforward algebra one 

gets U(\$+S$\) = C/( v /((^)+<5^) 2 +^2 +8<%) ~ 

U((<f> z )) + U'MJW, + (l7'«0,))/2<0,))(ty£ + H 2 ) + 
(U" ({(f) z )) /2)8<f) 2 . The term linear in 8<f) z is precisely 
compensated by —\/~ZH s 5(j) z coming from the external 
field term. The above expression for the quadratic terms 
in the expansion of U(\(j> + 6<j>\) should not come as 
a big surprise. Indeed, for the isotropic system the 
transversal staggered susceptibility in the presence of a 
finite staggered field H s is given by Xx = M S (H S )/H S , 

(s) 

while the longitudinal staggered susceptibility is x\\ — 
dM s (H s )/dH s . On the other hand, in the single mode 

approximation one has x« \H S ) = Zv/A^ (H s ). This ar- 
gument again fixes the longitudinal and transversal en- 
ergy gaps at 
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At this point it is important to see what assumptions 
the derivation of Eqs.(|7, ||) relies on. Indeed, we implic- 
itly assumed that the only parameter of the system which 
changes in the presence of staggered field are energy gaps 
Aj_ , and Am. We have disregarded the changes in Z and 
v with field. It is well known that the spin wave veloc- 
ity is not very sensitive to the parameters of the system 
and changes only slightly when, for instance, anisotropy 
is switched on Jl3| . This is rather natural result in field 
theoretical formulation of the problem since the velocity 
should not be renormalized at all except due to intrin- 
sic microscopic asymmetry of spatial and temporal direc- 
tions (spatial coordinate has underlying discrete lattice 
structure, while time is naturally continuous). The weak- 
ness of the change of Z with the field is due to the fact 
that in two dimensions the critical exponent rj, relating 
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Z to the correlation length £ as Z ~ is equal to zero. 
Connecting Z to its NLSM counterpart go in ([!]) one can 
approximately write £ ~ exp(— 2n/Z) or, alternatively, 
Z" 1 ^) = const - (2vr)- 1 ln£ = const + (27r) _1 In A. We 
see that a change in the gap A leads only to the loga- 
rithmic change in Z. Such corrections are disregarded 
at the precision of our calculations. In (^) the weak de- 
pendence of Z on £ is reflected in the omission of the 
nonlinear terms containing derivatives of <fi. It is a well 
know fact |ll| that for small values of r\ this is a valid 
approximation. We conclude with the derivation of an 
explicit expression for the staggered field dependence of 
longitudinal and transversal gaps for 5 = 1 Heisenberg 
chain. Plugging the numerical values for (74, 176: 58 1 an d 
Z in Eqs. (jj, U) we get 



Aj(M s 

A 2 
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= 1 + 1.56A/ 2 + 2.4M 4 + 3.27M S 6 (9) 



1 + 4.68A/ 2 + 12M S 4 + 22.9M S 6 (10) 



[19] J. Berges, N. Tetradis, and C. Wetterich, Phys. Rev. Lett. 
77, 873 (1996). 



FIG. 1. The staggered magnetization curve as de- 
duced from neutron powder diffraction data Q) on 
(Nd^Yi-^BaNiOs (x = 0.25,0.5,1). The effective stag- 
gered field is estimated through a mean-field analysis [[| . The 
solid line is a single-parameter fit with our theoretical result 
X M (0)H s = M s (l + 1.56M S 2 + 2.4M 4 + 3.27M 3 6 ) (see text). 



FIG. 2. Inelastic neutron scattering data for the relative 
increase in the energy gap of the Ni-chain magnon as a func- 
tion of induced staggered spin M 3 in Nd^Yi-^BaNiOs jsjj. 
Solid lines are defined by (9,10) The experimental data agree 
with the theoretical prediction for the transversal gap. 



As shown in Fig. 2, our prediction for the transversal gap 
is in excellent agreement with neutron scattering data on 
(Nd x Yi_ x )2BaNi05 family of compounds |ict| , with no 
adjustable parameters. 

The work at Brookhaven National Laboratory was sup- 
ported by the U.S. Department of Energy Division of 
Material Science, under contract DE-AC02-98CH10886. 
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